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We study the electronic and magnetic structures of the newly-discovered compound KCrsAss. The 
non-magnetic state has five Fermi surface sheets involving respectively three quasi-one-dimensional 
and two three-dimensional energy bands. However, the ground state is magnetic, exhibiting a novel 
interlayer antiferromagnetic order where the basic block-spin state of a unit Cr triangle retains a high 
spin magnitude. Moreover, its Fermi surface involves three one-dimensional sheets only, providing 
evidence for local moments in this compound due to the reduced dimensionality. By fitting a twisted 
spin tube model the magnetic frustrations caused by local moments are found to be relaxed, leading 
to gapless spin excitations. A frustration-induced transition to the disordered low block-spin state 
is expected upon increasing the intralayer exchange interaction. 


PACS numbers: 75.10.Pq, 75.50.Lk, 75.70.Tj 

Introduction. The interplay between superconductiv¬ 
ity and magnetism or other competing density-waves 
has been a major theme in modern condensed mat¬ 
ter physics, particularly in transition metal compounds 
with both strong electron correlations and reduced di¬ 
mensionality. Although the development of long-range 
orders at finite temperature is forbidden in strictly 
one-dimensional system with short-range interactions^’^, 
such restriction can be lifted for quasi-one-dimensional 
(QID) systems. Typical examples include the QID su¬ 
perconductors (TMTSF)2X^, M2Mo6Se6 (M=Tl,In)"^, 
Lio.gMoeOir^’®, etc. A ’’universal” phase diagram for the 
QID superconductors suggests that the superconductiv¬ 
ity is in close proximity to the density-waves 

The recently discovered QID superconductors 
AgCraAss ( A=K, Rb, Cs) provide a new setting for 
understanding the interplay of unconventional super¬ 
conductivity and density wave instabilities in reduced 
dimensions.The building blocks of these compounds 
consist of QID [CrAsJoo double-wall nano-tubes inter¬ 
calated by A+ cations. The first principle calculations 
using the density functional theory (DFT)^^’^^ predicted 
a three-dimensional ( 3 D) Fermi surface (FS) sheet and 
two QID FS sheets, all contributed mainly by the 
Cr- 3 (i electrons. The penetration depth measurement^^ 
has evidenced a nodal line in the pairing state below 
Tc- Meanwhile, the NMR experiments^ has revealed 
a power law behavior of the spin-lattice relaxation 
rate, manifesting the ID Luttinger liquid feature at 
the normal state. A couple of theoretical proposals 
have been devoted to the possible pairing symmetries 
of the superconducting state, concerning the respective 
relevance of the QID and 3 D bands^^’S'^^ while the 
formation of molecular orbital bands has been detailed 
in a twisted Hubbard tube.^® These studies suggest 
possible competing superconducting and density wave 
ordering tendencies in K2Cr3As3. In fact, supercon¬ 
ductivity associated with a magnetic quantum critical 
point has revealed in the bulk materials of chromium 


arsenides^®’^®. The material realization of density-wave 
instabilities competing with superconductivity in the 
QID chromium arsenide compounds should be very 
interesting and require more surveys. 

Motivated by this issue, we will study the elec¬ 
tronic and magnetic structures of a related compound, 
KCrsAss. This 133 -type of chromium arsenide has been 
recently synthesized experimentally^^. It consists of the 
QID [CrAsJoo structure similar to K2Cr3As3. The com¬ 
pound is metallic down to the lowest measured temper¬ 
ature, but no superconductivity has been observed. In¬ 
stead, the polycrystal sample displays a spin glassy fea¬ 
ture at low temperatures^^. In this paper, we find based 
on DFT calculations that the electronic and magnetic 
structures of KCr3As3 exhibit new properties in con¬ 
trast to the superconducting 233 -type compound. The 
most prominent features are the lack of three-dimensional 
Fermi surface sheets and emergence of the interlayer anti¬ 
ferromagnetic order. These numerical results indicate co¬ 
existence of local moments and itinerant electrons in the 
KCr3As3 compound due to the reduced dimensionality. 
Such exotic metallic magnetic phase is understood as due 
to the formation of high spin states in each Cr triangles 
via the orbital selective Mott transition. A frustration- 
induced quantum phase transition to the disordered low 
block-spin state with a spin gap is expected by tuning the 
intralayer exchange to the antiferromagnetic side. Our 
study thus reveals a possible interplay between uncon¬ 
ventional superconductivity and magnetic order in this 
class of materials. 

Calculation Method and Details. All reported re¬ 
sults are obtained using the DFT as implemented in the 
Vienna Abinit Simulation Package (VASP)^^’^®, where 
the valence electron-ion interactions are modeled us¬ 
ing projected augmented wave (PAW) method^®’^^ and 
the exchange-correlation effects are approximated with 
Perdew-Burke-Enzerhoff flavor of general gradient ap¬ 
proximation (PBE).^® To ensure the convergence of the 
calculations, a 540 eV energy cutoff to the plane wave 
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basis and a 6 x 6 x 12 F-centered K-mesh is employed; 
whereas a 12 x 12 x 24 F-centered K-mesh and tetrahe¬ 
dron method are employed to perform the DOS calcula¬ 
tions. 




FIG. 1: Crystal structure of KCrsAsa. (a) Experimentally 
determined geometry with space group Pda/m (No. 176); (b) 
Hypothetic geometry closer to K 2 Cr 3 As 3 with space group 
P6m2 (No. 187). In both panels, the red, blue and green balls 
represents As, Cr and K, respectively; larger and smaller balls 
indicate atoms at z=0.75 (X^) and z=0.25 (X^^), respectively 
(X=K, Cr, and As). 

Experimentally, the crystal structure of KCrsAsa 
shown in FIG. 1(a), was determined to be centrosym- 
metric with space group PBa/m (No. 176). Unlike the 
K 2 Cr 3 As 3 structure, the Cr^-Cr^ bond lengths are iden¬ 
tical to the Cr^^-Cr^^ bond lengths in KCr 3 As 3 . Apart 
from this apparent difference, the [CrAsJg sub-nanotubes 
in KCr 3 As 3 are rotated ~ 10 ° with respect to their cen¬ 
tral axis compared to the K 2 Cr 3 As 3 structure. In or¬ 
der to make connections with the 233-compound, we 
have also performed calculations on a hypothetical crys¬ 
tal structure of KCr 3 As 3 which is initially with space 
group PGs/m (No. 176) but without the rotation of the 
[CrAsJe sub-nanotubes (FIG. 1(b)). Interestingly, this 
structure spontaneously turns into P6m2 (No. 187) af¬ 
ter optimization. So the loss of symmetric center is not 
due to the distribution of K atoms between z = 0.25 and 
z = 0.75 planes, but rather relates to the rotation of the 
[GrAsJe sub-nanotubes. 

Table I lists the total energies of KGr 3 As 3 with differ¬ 
ent magnetic configurations. Unlike the K 2 Cr 3 As 3 com¬ 
pound, where multiple magnetic states are energetically 
almost degenerate in DFT calculations, the inter-layer 
antiferromagnetic (lAF) configuration, a kind of block 
spin magnetic order, is apparently energetically much 
lower and is stable against structural optimization. Also, 
it is worthy noting that both the non-collinear phases 
(lOP and NCL) are energetically much higher than the 
collinear phases (FM and lAF). Upon optimization, the 
collinear phases are energetically stable while the non- 
collinear ones are not. For the lowest energy state lAF, 
the magnetic moment on each Cr is calculated to be 
~0.77 /is, in good agreement with the experimental ob¬ 
servation of 0.68 

In order to understand the formation of magnetic 
moments, we first examine the electronic structure of 
KCr 3 As 3 in the non-magnetic phase. As shown in FIG. 
2 (a), the band structure of KCr 3 As 3 near the Fermi 
level is quite similar to that of K 2 Cr 3 As 3 , including 
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FIG. 2: (a) Band structure and (b) DOS of non-magnetic 
KCr 3 As 3 . In panel (a), the size of the red/blue/green cir¬ 
cles are proportional to the orbital weight of the Cr-3d22/- 
?)d^2_y2/-^dxy, respectively. The dashed horizontal line in¬ 
dicates the Fermi level corresponding to doping the system 
with 1 electron/formula, i.e. K 2 Cr 3 As 3 . 


the orbital characters, but the degeneracies at F and 
K are completely different. Indeed, the degeneracies at 
K of KCr 3 As 3 and the double degeneracy of the bands 
around Ep-Q.b eV at ^ 2 = 0 .5 (A-H-L-A) are direct con¬ 
sequences of the P 63 /m symmetry which is broken in 
K 2 Cr 3 As 3 . Furthermore, the reduction of 50% potas¬ 
sium atoms causes hole doping to the system, resulting 
in a shift of the Fermi level. As a result, the number of 
bands cross the Fermi level becomes five, among which 
three are quasi-lD and two are 3D. The DOS and the 
orbital contributions are shown in FIG. 2(b). The quasi- 
ID signature van Hove singularities are apparent around 
-1.6 eV, -0.5 eV, -0.4 eV, 0.5 eV relative to Ep. The DOS 
around Ep is dominated by the Gr-3d orbitals, in partic¬ 
ular, Cr-3da;y and 3 da; 2 _j ,2 orbitals. The contribution of 
Gr- 3 d ^2 orbital is small, but non-negligible since the a- 
band is almost exclusively contributed by Cr- 3 d 22 (FIG. 
2(a)). The Fermi surface (FS) of KCr 3 As 3 constitutes 
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TABLE I: Total energies of KCraAss with different magnetic confignrations. Column Expt, Opt and P6m2 list the results with 
experimental lattice constants and coordinates, the optimized lattice constants and coordinates, and the hypothetical structure 
after optimization, respectively. Please be noted that the hypothetical structure was initially centrosymmetric Pha/m (No. 
176), but spontaneously turned into P6m2 (No. 187) after optimization. 




EM 1 lAF |IOP|NCL| NM 

BIAF= 1UUD'' 1 FM= | lAF'' | NM'^ 

Exp 

a (c) 

AE (meV) 
mcr (ps) 

-41.3 (-43.9) 
0.36 

9.0886 (4. 
-121.7 (-119.4) 
0.77 

1801) 

-22.5 -16.0 
0.51 0.31 

0.0 

0.00 

9.0886 (4.1801) 

-20.7 -23.4 -50.7 -129.7 0.0 
0.44 0.40 0.37 0.78 0.00 

Opt 

a (c) 

AE (meV) 
mcr (ps) 

9.3442 (4.0862) 
-0.8 ^6.0) 
0.14 

9.3419 (4.0974) 
-20.4 

0.29 

to NM 
to NM 
to NM 

9.3447 (4.0826) 
0.0 

0.0 


P6m2 

a (c) 

AE (meV/Cr) 
mcr (ps) 

9.7091 (4.0848) 
-1.6 (-4.1) 

0.1 

9.7067 (4.0904) 
-7.0 

0.24(-0.29) 

to NM 
to NM 
to NM 

9.7109 (4.0826) 
0.0 

0.00 



of three quasi-ID sheets and two 3D sheets as shown in 
the Supplemental Material (SM). If the system is doped 
with one electron per formula (i.e. K 2 Cr 3 As 3 ), the FS 
sheets will be reduced to two quasi-lD and one 3D, as in 
the original K 2 Cr 3 As 3 compound®. Of course, the shape 
of the 3D FS sheet will be different from the one in the 
original K 2 Cr 3 As 3 due to the local structure difference. 

We now show the electronic structure of the lowest in¬ 
terlayer antiferromagnetic phase (lAF) in Figure 3. The 
orbital characteristics of bands near Ep remains the same 
as those in the NM state. The four Cr-3d3;2_j,2/Cr-3da;y 
bands near Ep are separated by ^0.4 eV in the lAF 
phase, presumably due to formation of local moment on 
Cr; whereas these bands are very close to each other in 
the NM phase. As a result, two of these bands, who 
formed the 3D FS sheets in the NM phase, become fully 
occupied, leaving us only three quasi-lD sheets in the 
lAF phase. It is worthy noting that the band around 
Ep-0.03 eV are very flat, yielding strong quasi-lD van 
Hove singularity in the DOS plot. These electrons may 
give rise to activation-like behavior in resistivity^^. 

Local Moments and Magnetic Frustration. The above 
DFT calculations suggest the existence of local moments 
in CrsAss, given the facts that the calculated magnetic 
moment 0.77 ^b is close to the saturated value of 
S' = 1/2 isolated local moment and quantum frustra¬ 
tions are not considered so far. The local moment phase 
is also hinted in the non-magnetic band structure FIG. 2 
(a), where two upper degenerate bands along the F —>■ A 
direction intersect the Fermi energy at the Fermi point 
kp^ O.SGtt. This value is very close to 7r/2, indicating 
the nearly half-filling case of the corresponding molecular 
orbital band. In ID, any repulsive Coulomb interaction 
will induce the 4:kp^ Umkllamp scattering at half-filling, 
resulting in a finite charge gap. Therefore, a Mott transi¬ 
tion in this band is naturally expected in the correspond¬ 
ing molecular orbital band^®. In fact, compared with the 
band structure of K 2 Cr 3 As 3 , both the band narrowing 
effect and the electron correlations are enhanced due to 
the reduced dimensionality in KCr 3 As 3 , providing more 
evidence for a QID Mott phase via orbital selective Mott 
transitions^®’^^ as in iron pnictides®®^®®. 
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FIG. 3: (a) Band structure and (b) DOS of KCrsAsa at the 
lAF state. Only spin up channel is shown for clarity pur¬ 
poses. In panel (a), the size of the red/blue/green circles are 
proportional to the orbital weight of the Cr-3d,j2/-3d„,2_j,2/- 
3dxy, respectively. The dashed horizontal line indicates the 
Fermi level corresponding to doping the system with 1 elec¬ 
tron/formula, i.e. K 2 Cr 3 As 3 . 
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We are then led to consider a local moment at each 
Cr site described by a spin operator, Sr, with the mag¬ 
nitude S = 1/2. Assuming that the coupling between 
itinerant electrons and local moments is relatively weak, 
we consider the Hamiltonian of these localized spins as 

^spin = Jl ^n,I ' ^n,I' + Jl ^n,I ' ^n,I' (1) 

n,I.r nJJ’ 

+ J2 E Sji,/ • // -|- J2 ^ ^ 

n,lj' n,I,r 

Where, the first two terms indicate the intra-layer 
nearest-neighbor couplings in the two different triangles, 
while the last two terms the inter-layer nearest-neighbor 
couplings between the two triangles within one unit cell 
and across two unit cells, respectively. 

The spin model is an extended variant of twisted 
spin tubes previously proposed for the three leg systems 
[(CuCl 2 tachH) 3 ]Cl 2 and CsCrF 4 .^^“^‘* The magnetic frus¬ 
trations will arise not only due to the intralayer spin 
exchanges but also the interlayer or inter-unit-cell spin 
exchanges if they are all antiferromagnetic. Concerning 
the present KCraAss compound, we can assume Ji = Ji, 
J 2 = J 2 . Their values are estimated by fitting the cor¬ 
responding energies in different magnetic structures with 
the moment per Cr atom being fixed at 0.77 hb- The 
best fitting suggests Ji = —9.9 meV and J 2 = 18.5 
meV. It is worthy noting here that we have also deter¬ 
mined the strength of the inter-tube interactions Jt for 
the lAF state to be < 0.1 meV by comparing the DFT to¬ 
tal energies of inter-tube-FM and inter-tube-AFM states. 
Thus the KCrsAsa system is a truly remakable exam¬ 
ple of nearly one-dimensional spin chain, explaining why 
the lAF long range order was not observed in experi¬ 
ments. The ferromagnetic Ji implies the suppression of 
the magnetic frustration in this 133-compound. In this 
case, the twisted tube is basically a rhombus lattice with¬ 
out frustration due to the bipartite nature of the lattice. 
A large block-spin S'a = 3/2 state is then expected for 
each Cr triangle, involving the four polarization compo¬ 
nents: If) = I ttt), I - I) = I I5) = tu 
) + I t/t) + I itt)]) I ~ 5) = "^[l + I iti) + I iit)]- 

By introducing a spin-3/2 operator T acting on these 
Sa = 3/2 states, we obtain the effective Hamiltonian 

'Hsa=3/2 = E[^^ "Tn ' "Tn+l + K2{Tn- 'Tn+1 )"] (2) 

n 

with Ki = 2J2/3-bJ|/(18|Ji|) and A :2 = - J|/(54| Ji|)35. 
The low temperature property is dominated by the Ki- 
term. With quantum fluctuations the interlayer antifer¬ 
romagnetic order is only quasi-long-ranged, in the sense 
that the spin-spin correlation function exhibits a power- 
law behavior and the spin excitation spectrum is gapless. 

This magnetic ground state remains stable in a wider 
parameter regime when Ji moves toward the antiferro¬ 
magnetic side. However, when Ji is much larger than 
J 2 , the ground state will be changed due to strong frus¬ 
trations. In the limiting case, each triangle can be 


solved independently, giving rise to the doubly degen¬ 
erated ground states involving | ts) = ■;^[| tti) +w| 

) + I ti) = ;^[| m) + tit) + ttt)] 

(and the states with | f) <^ | J,)). Noted that these lower 
block-spin Sa = 1/2 states are chiral in nature, associ¬ 
ated with the C 3 rotation w = To the lowest order 

in J 2 , the effective spin model describing these Sa = 1/2 
spin excitations is given by^®’^^ 

^Sa=i /2 = E ’ Tk+i[1 -I- 2(a;^T/~r„^3 -I- h.c.)], 

n 

with T being a spin-1/2 operator acting on the two dif¬ 
ferent spin polarization components cr =t, t; while 
on the chirality degrees of freedom, namely, r“'"|crL) = 
0,T“|(Ti?) = |L),T+|cri?) = \L),T~\aR) = 0. According 
to the early numerical studies, this model has a disor¬ 
dered ground state, i.e., the dimerized ground state ex¬ 
hibiting gaps in both spin and chirality excitations. 

It is well-known from the celebrated Lieb-Schultz- 
Mattis theorem^® that half-spin chains have either a de¬ 
generated ground state or gapless spin spectrum. Inter¬ 
estingly, with increasing ratio a = J 2 IJ 1 , both situations 
could be realized in the present material. The quantum 
phase transition from gapless to gapful phases takes place 
around Uc ~ 1.22®®’®9. 

Summary and discussions. Our DFT calculations 
show that the ground state of the KCraAsa compound is 
metallic and extremely QID dimensional compared with 
its superconducting 233-counterpart. Moreover, it ex¬ 
hibits an interlayer antiferromagnetic order with a large 
block-spin 3 x 0.77 ^b ~ 2.3/iB in each Cr triangles. By 
fitting the twisted J 1 -J 2 spin tube model, we find that 
the basic magnetic property can be actually captured 
by the quantum antiferromagnetic Heisenberg chain with 
block-spin Sa =3/2 defined in Eq.(2). The absence of 
superconductivity is likely due to the strong competition 
tendency of the quasi antiferromagnetic order in this 133- 
compound. 

We should remark that due to the reduced dimen¬ 
sionality, a static long-range magnetic order is unstable 
against to enhanced quantum or thermal fluctuations, 
and strongly influenced by the disorder effect. So it is 
not very strange that the spin glass behavior, instead 
of the static magnetic order, is observed in the polycrys¬ 
talline samples^^. We also remark here that although the 
ground state itself is gapless, the lower spin Sa = 1/2 
sector may still play some role at finite temperatures, ex¬ 
hibiting a crossover-to-spin-gap behavior as in the spin 
tube material [(CuCl 2 tachH) 3 ]Cl 2 ^°. 

An intriguing situation with a large intralayer anti- 
ferromanetic coupling Ji is expected, hopefully by ap¬ 
plying the planar pressure while keeping the symmetry 
of charge distributions among the two conjugated trian¬ 
gles. Such situation should be contrasted to the 233- 
compound, where the more K concentration induces the 
non-centrosymmetric structural distortion and asymmet¬ 
ric charge distributions in the two conjugated triangles. 
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Because superconductivity may be also induced from the 
spin-gaped phase upon doping'^^ and given the recent 
pressure experiments on the K 2 Cr 3 As 3 ^^’^^, we expect 
that the disordered phase with both spin and chirality 
gaps can be realized in further coming pressure experi¬ 
ments on KCr 3 As 3 and the related compounds. 
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SUPPLEMENTARY MATERIAL 

Reduced Dimensionality and Magnetic Frustration 
in KCrsAss 


FIG. SI: Fermi surfaces of (a) prestine KCraAsa and (b) 
electron-doped KCrsAsa. 


by: Chao Cao, Hao Jiang, Xiao-Yong Feng, and Jian- 
hui Dai 



(a)Prestine KCraAsa 




The Fermi surfaces are reconstructed by fitting the 
first-principles band structure to a 10-orbital tight- 
binding Hamiltonian using the maximally localized Wan- 
nier functions (MLWFs). Figure Sl(a) shows the Fermi 
surface sheets of KCraAsa. Three quasi-lD sheets and 
two 3D sheets constitute the Fermi surface of prestine 
KCraAsa as expected. If the system is doped with le per 
formula (i.e. K 2 Cr 3 As 3 ), the FS sheets will be reduced 
to two quasi-lD and one 3D, as in the original K 2 Cr 3 As 3 
compound^^. However, the shape of the 3D FS sheet 
will be completely different from the one in the original 
K 2 Cr 3 As 3 due to the local structure differences. 



















































































